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ABSTRACT 

We use Monte Carlo realizations of halo formation histories and a spherical 
accretion model to calculate the expected scatter in the velocity dispersions 
of galactic halos of a given mass due to differences in their formation times. 
Assuming that the rotational velocity of a spiral galaxy is determined by the 
velocity dispersion of its halo and that its luminosity is related to its total 
baryonic mass, this scatter translates to a minimum intrinsic scatter in the Tully- 
Fisher relation. For popular cosmological models we find that the scatter due to 
variations in formation histories is by itself greater than allowed by observations. 
Unless halos of spiral galaxies formed at high redshift (z ^ 1) and did not later 
accrete any significant amount of mass, the Tully-Fisher relation is not likely to 
be the direct result of cosmological initial conditions but rather a consequence of 
a subsequent feedback process. 
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1. Introduction 

One of the most surprising observational results in extragalactic astronomy is that the 
luminosity L and rotational velocity v c of a spiral galaxy obey a tight power-law relation, 
L oc v", with a ranging between a value of 2 in the blue and a value of 4 in the infrared 
(Willick & Strauss 1995, and references therein). Because v c can be directly measured while 
the inferred value for L depends on distance, this so-called Tully-Fisher relation (Tully & 
Fisher 1977; Aaronson, Huchra, & Mould 1979) has been widely used as a tool to measure 
distances. By now, various catalogs have been compiled with an interest in calibrating and 
applying this relation. Differences in observational technique, such as bandpass, the method 
for observing the rotational velocity, and sample selection, affect the power-law slope and 
degree of scatter in the relation, with I band observations providing the tightest fit and 
a 2.7 (Mathewson et al. 1992). In the largest samples the scatter in the relation is around 
0.4 mag over two decades in luminosity (Mathewson & Ford 1994; Willick et al. 1995a, b,c), 
while in other samples, scatter as low as 0.1 mag has been reported (Bernstein et al. 1994). 
Because of the correspondingly low distance errors, the Tully-Fisher relation has been widely 
used in peculiar velocity studies (Strauss & Willick 1995 and references therein). 

Despite the popularity of the relation, the reason behind its existence is still unknown. 
It is particularly surprising that the complex processes of gravitational collapse and star 
formation do not act to blur this tight relation. In this paper we address the most basic 
question in this context, namely: Could the low scatter observed in the Tully-Fisher relation 
be the direct result of cosmological initial conditions? We therefore wish to relate the two 
quantities in the relation, the circular velocity and the luminosity, to properties of the dark 
matter halos in which the spiral galaxies form. The former depends on the shape and 
depth of the halo potential-well, while the latter depends on the total baryonic mass. If 
the cosmological scatter in these properties is greater than allowed by observations, then 
the Tully-Fisher relation would imply a subsequent feedback process that regularizes gas 
dynamics and star formation in galaxies according to the depth of their potential wells. 

In this paper we study the spread in the potential depth of an ensemble of galaxy halos 
of a given mass due to the differences in their formation times. We construct Monte Carlo 
realizations of the formation histories of halos as provided by the excursion set formalism 
(Bond et al. 1991, hereafter BCEK; Lacey & Cole 1993, 1994, hereafter LCa,b). The history 
of a given halo is used to estimate its current binding energy and velocity dispersion under the 
most regularizing assumption that the halo relaxes to an isothermal configuration through 
spherical accretion. This method and its minimum scatter estimates are described in §2 for 
different halo masses and cosmological models. In §3 we relate these results to the observed 
scatter in the Tully-Fisher relation. We conclude that the implied minimum scatter due to 
the variations in halo formation histories is too high to be reconciled with observations. 
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2. Method and Results 

The universality of flat rotation curves and the Tully-Fisher relation imply that spiral 
galaxies obey a regular behaviour. In order to show that this regularity cannot naturally 
emerge from a random set of initial conditions, it is sufficient to consider the most optimistic 
model in which all galaxy halos form as isothermal spheres. Variations on this assumption 
would only tend to smear further any regularity we find. Collisionless N-body simulations and 
analytic results indicate that halos do indeed tend toward distributions resembling isothermal 
spheres (Dubinski & Carlberg 1991; Katz & Gunn 1991; Spergel & Hernquist 1992). 

The potential depth of an isothermal halo is simply characterized by its one-dimensional 
velocity dispersion a. If all, or some fixed fraction, of the baryons in the halo cool to the 
galactic disk, then a perfect relation between the velocity dispersion a and the mass of 
the halo M would produce virtually no scatter in the Tully-Fisher relation. This follows 
from the most optimistic assumption that the circular velocity of the disk at large radius 
is simply proportional to the velocity dispersion of the halo, neglecting any complications 
due to triaxiality of the halo, ellipticity of the galactic disk (Franx & de Zeeuw 1992), or 
variations in the rotation curve due to self-gravity of the baryonic disk. 

However, as illustrated by the spherical collapse model (Gunn & Gott 1972; White 
1995), a is not simply a function of M. Because a 2 characterizes the kinetic energy of a 
virialized system, we can fix a by requiring that energy be conserved through the virialization 
process. For the collapse of a spherical shell without a cosmological constant (A = 0), the 
shell energy per unit mass is dE / dM = —GM/R ta} where M is the mass enclosed within the 
shell and R ta is the radius of the shell at turn-around. For A = cosmologies, the spherical 
shell solution gives 

Rta oc M^tlH (1) 

where t co a is the time of collapse. We may then integrate over all shells to find the total 
energy E, which yields the velocity dispersion a 2 oc E/M. If an object collapses at a single 
time, a oc M 1 ^H C ^ 3 . The constant of proportionality is unimportant for this work. 

During the hierarchical process of structure formation, the formation of an object is not 
characterized by a single redshift of collapse. Rather, the object is assembled over time by 
a combination of mergers and smooth accretion. Thus, one would expect variations in these 
formation histories to produce a scatter in the velocity dispersion for an ensemble of objects 
of a given mass. To approach this problem analytically, we apply the excursion set formalism 
(BCEK) that has been used to predict halo formation histories (LCa,b). Given a collapsed 
object of mass M 2 at some redshift z 2 this formalism predicts the probability distribution 
for the object to have a mass Mi at some higher redshift z\. We denote by (Tm{M) the 
rms fluctuation amplitude 6M/M (using a k-sp&ce top-hat window function) as a function 
of mass M and define 6 c (z) to be the overdensity required for the collapse of a spherical 
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perturbation to occur at redshift z\ both these quantities are measured at the present epoch 
as extrapolated in linear theory. Then the value of M\ is found implicitly from the relation 



where x is a Gaussian-distributed random number with zero mean and unit variance (see 
eq. 2.16 of LCb and change variables to x 2 = [6 c i — 6 c2 ] 2 /[erf — cr|]). Using this equation, we 
start from a given mass at the present epoch and step backwards in time to find its mass 
at a sequence of higher redshifts. We continue this procedure until Mi is less than 10% of 
the present-day mass. Note that because we start from the present-day mass of the object, 
we are restricting our sample to isolated galaxies. Because isolated spiral galaxies do satisfy 
the Tully-Fisher relation, this is an acceptable restriction. 

Next, we need to convert a given mass history into a binding energy. We do so by 
assuming that the mass gained in an infinitesimal redshift interval was in fact accreted in 
the most regular fashion, i.e. in spherical symmetry and with a constant dM/dz. Any 
variation on these simplifying assumptions would introduce additional degrees of freedom, 
such as the orbital parameters or the internal structure of the merger components. For 
spherical accretion in an arbitrary cosmology, the infinitesimal energy added to the system 
is 



where fi is the present fraction of the critical density in non-relativistic matter, Oa is the 
present fraction of the critical density in the cosmological constant, = 1 — fi — at 
the present epoch, M is the mass interior to the added shell, dM is the shell mass, and 



The last factor in equation 3 is the generalization of the factor of t coll in equation 1. For 
finite Az and AM, we integrate over AM assuming dM/dz = AM/ Az to find AE. We 
then sum over all the redshift steps to find the total energy, which we set proportional to 
the square of the velocity dispersion. 

As discussed in LCa, the above merger histories allow the possibility that the mass 
accreted in a given time step (AM) may be greater than the mass in the progenitor (Mi). 
This is in conflict with our goal to track the mass of the largest progenitor of an object. In the 
context of equation 3, it means that we are allowing most of the object to accrete spherically 
onto a small central seed in a short redshift interval. As this process seems unrealistic, we 
reverse the roles of the two pieces and consider a small spherical accretion onto a larger mass 
central piece. In each redshift step Az, if the mass Mi found from equation 2 is less than half 






(3) 




(4) 



- 5 - 



Fig. 1. — Fractional scatter in the velocity dispersion a as a function of the present-day 
overdensity 6M/M = (Tm(M ) on the mass scale of the final object M . The different curves 
correspond to scale-free power spectra, P(k) oc k n , with different values of n, in an fi = I 
cosmology. The "No Mergers" case requires that the object never accrete more than 20% of 
its final mass in any given time step. 

of Af 2 , we flip the roles of the pieces and instead take the mass to be M 2 — M\. LCa use a 
slightly different scheme to remedy this problem when constructing their Monte Carlo merger 
history algorithm. Surprisingly, when they compare the formalism to N-body simulations 
(LCb), the statistics of the time at which a present-day object has acquired half its mass 
better match the unconstrained formalism rather than the Monte Carlo formalism, despite 
the fact that following the more massive of the two merging pieces seems more appropriate. 
In this work, we use the constrained histories (i.e. tracking the more massive piece), since 
they produce about 15% less scatter in the M-a relation than the unconstrained histories. 

We may now integrate over many generated mass histories to find the variation in 
the velocity dispersion on a given mass scale. We first consider scale-free power-spectra 
P(k) oc k n , for which om oc Af ~( n+3 )/ 6 , and find Acr/cr as a function of <tm(Mo) and n, where 
Mo is the present-day mass and Act is the standard deviation of the distribution of velocity 
dispersions a. The results are shown in Figure I and compared to the scatter in luminosity 
that would result if we enforced a Tully-Fisher relation of the type L oc a 3 (dotted horizontal 
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Fig. 2. — Fractional scatter in the velocity dispersion a in four different CDM models. The 
flat, low-fio model is the only one with a non-zero cosmological constant. The parameter 
T = Cl h is used to fit the power spectrum (cf. Efstathiou et al. 1992). 

lines). The figure shows that the scatter in a drops strongly with decreasing (Tm(M ) and 
with more negative n. Both these changes tend to bring the collapse to lower redshift so 
that the ensemble of objects is described by a narrower range of collapse redshifts. An open 
universe (O < 1, A = 0) results in a larger scatter at a given (Tm(M ) and n. A flat universe 
with A / produces the same scatter as its flat A = counterpart, because if 0# = then 
both equations 2 and 3 involve only 8 C and functions of mass, with A entering only in the 
constant of proportionality of the energy and therefore dropping out of A<t/<t. 

We next consider various cold dark matter (CDM) models, as parameterized by the 
constant Y (Efstathiou et al. 1992). We normalize the power-spectrum to COBE using <r 8 , 
the rms amplitude of fluctuations in spheres of radius 8h _1 Mpc (Gorski et al. 1995a, b; 
Stompor et al. 1995). Figure 2 presents the resulting scatter as a function of mass scale. We 
find that = 1, COBE-normalized CDM (solid line) produces too much scatter; this reflects 
the tendency of this model to have too much power on small-scales. Under-normalized 
= 1 CDM (dotted line) produces less scatter, while open and A-dominated universes give 
yet lower values. The lower value of <r 8 is the major factor in reducing the spread. However 
lower values of <r 8 than the ones shown in the figure are not allowed by the abundance of 
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clusters of galaxies (White et al. 1993). 

Finally, we consider the actual M-a relation rather than its expected scatter. If all of 
the objects had the same formation history regardless of mass scale, then E oc Ma 2 oc Af 5 / 3 , 
would yield M oc a 3 . Instead we find a weak tendency for smaller objects (with higher 
<7m[Mo]) to collapse at higher redshifts and therefore to have more binding energy. We 
therefore find M oc a a , with a 3.1 — 3.2. Although a weak dependence of M/L on M may 
lower the expected exponent by a few tenths, this relation is in reasonable agreement with 
the observed Tully-Fisher relation in the R and I bands (Mathewson et al. 1992) as long as 
its scatter is ignored. 

3. Discussion 

We next examine the implications of the predicted scatter in M-a for the Tully-Fisher 
relation. In attempting to rule out a primordial origin for this relation, we adopt the most 
optimistic approach for connecting M-a to L-v c . We assume that the circular velocity of 
the baryonic disk is strictly determined by the velocity dispersion of its spherical isothermal 
halo and that the overall luminosity of the galaxy is tightly related to the mass of its halo. 
The latter relation might exist if all the baryons in the halo cooled onto the central disk and 
were turned into stars with a universal mass-to-light relation. We then treat the scatter in 
a on a given mass scale as a minimal source of intrinsic scatter in the Tully-Fisher relation. 

How much scatter is allowed by the observations? The full observed scatter in the largest 
samples is 0.4 ± 0.02 magnitudes (Willick 1995a, b,c). Many different effects contribute to 
this error budget. Observational errors are usually estimated to be 0.15-0.30 mag (Strauss 
& Willick 1995), with a significant contribution coming from the measurement of the HI line 
widths. The ellipticity of the halo potential leads to error in measuring the circular velocity 
and estimating the inclination of the baryonic disk (Franx & de Zeeuw 1992); these effects 
have been measured in a sample of spirals and estimated to contribute a scatter of about 
0.15 mag (Rix & Zaritsky 1995). The above effects alone reduce the remaining error budget 
to <^ 0.3 mag, which corresponds to about 10% scatter in the circular velocity. It seems 
likely that variations in the mass function of stars would lead to additional scatter, leaving 
less of the error budget available for intrinsic variations in a. In addition, other studies of 
the Tully-Fisher relation quote lower values for the observed scatter (Pierce & Tully 1988; 
Schommer et al. 1993; Bernstein et al. 1994). In summary, Aa / a = 10% provides an upper 
limit to the intrinsic scatter in a at fixed M that is allowed by observations. 

Naively, if each halo had a single formation redshift zj, then the above 10% limit would 
require that all halos of a given mass form within an unusually narrow redshift interval, 
Azj/(1 + zj) <^ 20%. Indeed, Figure 2 indicates that reasonable cosmologies are at best 
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marginally consistent with the above 10% limit. Smaller masses are predicted to have 
a larger scatter, matching a trend in the data (Federspiel, Sandage, & Tammann 1994), 
although various observational effects could also introduce larger scatter in dimmer galaxies. 

In calculating the minimum Tully-Fisher scatter due to differences in the formation 
histories of dark halos, we have introduced two simplifying assumptions. First, we assumed 
spherical accretion in relating the mass history of an object to its final binding energy E (cf. 
Eq. 3). Non-sphericity would add new degrees of freedom to our analysis, such as the orbital 
parameters or the internal structure of the merger components. These additional parameters 
could increase the resulting M-a scatter. Second, we have assumed that the mass and stellar 
luminosity of the gas that cools to form the galactic disk is a strict function of the mass of 
the halo. The ratio between the disk mass and the total mass of spiral galaxies is ^ 5—10% 
(Zaritzky & White 1994; Kochanek 1995). As this fraction is somewhat lower than the 
baryon fraction observed in clusters of galaxies (David, Jones, & Forman 1995), it is possible 
that not all of the baryons originally associated with a halo are currently in its stellar disk. 
If the fraction of baryons included in the disk varies from one galaxy to another, then the 
total luminosity of the galaxy will not be tightly constrained by the mass of its halo and 
the expected Tully-Fisher scatter would increase. As for the fate of the remaining baryons, 
they may either be expelled from the galaxy by supernova-driven winds (e.g. Spitzer 1956; 
Corbelli & Salpeter 1988) or remain in the halo in the form of hot pressure-supported gas or 
as dark compact objects. X-ray observations exclude the existence of extended x-ray halos 
or cooling flows in most spiral galaxies (Fabbiano 1989). Searches for microlensing events 
toward the LMC indicate that a noticeable fraction of the halo mass may consist of baryonic 
objects (Alcock et al. 1995; Gates, Gyuk, & Turner 1994). If these objects formed in the 
intergalactic medium at very high redshift, then they would be essentially indistinguishable 
from other types of cold dark matter. However, if they formed within proto-galaxies, then 
their mass fraction may vary from galaxy to galaxy. Similarly, if supernovae are effective at 
expelling gas out of the disk, then the amount of gas expelled would likely vary from one 
galaxy to another due to differences in their star formation history. Any such variations 
would contribute additional scatter to the relation. 

The Tully-Fisher relation applies to a particular morphological class of galaxies, and so 
one may argue that galaxies which have undergone major mergers since the formation of 
their stellar disks would not appear as thin spiral galaxies today (Toth & Ostriker 1992). 
This morphological restriction on merger histories ought to yield a smaller value of Acr/cr, 
but because the epoch of star formation is unknown, it cannot be incorporated in a simple 
way. Nevertheless, in order to examine the significance of this restriction, we have modified 
our algorithm so that any object that has undergone a merger that added more than 20% 
of its present-day mass is rejected. This constraint rejected 70% of the systems and reduced 
the scatter in v c by a factor of 20% (cf. the lowest curve in Fig. 1). Setting the rejection 
threshold at 30% in mass rejected 30% of the systems and reduced the scatter by 10%. While 
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this morphological constraint goes in the right direction, it leads to a small overall effect. 

Another way to implement the above morphological restriction is to argue that spiral 
galaxies already exist in their current form at some high redshift, zj. If accretion does 
not affect the rotation curve or luminosity of these galaxies between this early redshift and 
the present time, then the effective value of (Tm{Mq) would be lowered by the ratio of the 
growth factors at zj and today, e.g. by (1 + £/) in an = 1 universe. Note that the fact that 
structure formation cuts off at late times in an open universe is not sufficient for this purpose. 
As Figures 1 and 2 indicate, lowering (Tm{Mq) reduces the scatter in the M-a relation. At 
first glance, it seems unlikely that both the baryonic mass and the velocity dispersion of the 
halo would remain unchanged as the growth factor increases by a significant factor and the 
corresponding accretion occurs. However, detailed consideration of these conditions raises 
difficult questions concerning the maximum mass of galaxies (Thoul & Weinberg 1995), 
the accretion of satellite galaxies (Toth & Ostriker 1992; Navarro, Frenk, & White 1994), 
and the fate of dark matter cores as they fall into larger objects (Moore, Katz, & Lake 
1995). It is important to note that galaxy-sized halos in the universe at z = do not 
correspond to rare density peaks. Since gm ~ 3, they correspond to regions with a peak 
height v = S c (0)/a(M) ~ 0.5. However, galaxy formation arguments have often associated 
galaxies with v ~ 2 peaks. The excursion set formalism predicts that if galaxies form out of 
2cr peaks when the growth factor is one-third of its present value (i.e. z = 2 in an fi = 1 
universe or z = 4.3 in an open fi = 0.3 universe), then 90% (50%) of these galaxies are 
incorporated today within systems at least 4.4 (37) times as massive. Here we assume an 
n = — 1 power spectrum; more negative values of n result in even larger systems. Hence, 
although galaxy formation at high redshift reduces the M-a scatter, one must explain why 
the systems are unchanged by their infall into larger objects. 

Previous estimates of the M-a scatter caused by the properties of Gaussian random 
fields (Cole & Kaiser 1989) have assumed a Gaussian distribution of peak height v and then 
used a oc z/ 1 / 2 to derive the scatter in a. However, in the Press-Schechter formalism, objects 
of a particular mass at z = correspond to a single value of v. Peaks with larger v formed 
earlier and are now part of more massive objects. Since the Tully-Fisher relation is measured 
at z = 0, we cannot assign galaxies to a distribution of v. Instead, we must recognize that 
objects of a particular mass today were assembled through a variety of merger histories, 
leading to differences in their binding energies. It is interesting to note that these differences 
result in a smaller Tully-Fisher scatter (ss 0.4 mag) than that inferred from arguments based 
on the distribution of v (^1 mag). 

A different approach to the Tully-Fisher relation relies on the notion that the cooling of 
the gas is incomplete (Kauffmann, White, & Guiderdoni 1993, hereafter KWG). Here, one 
assumes that after collapse the gas acquires an isothermal density profile similar to that of 
the dark matter. If one then considers the cooling time as a function of radius and includes in 



- 10 - 



the galactic disk only the mass ol the gas that has had time to cool, Af co id, then the reservoir 
of cold gas depends on the velocity dispersion and not on the mass of the halo, since the 
former sets the density profile of the gas. Thus, it is possible to get a tight a — Af co id relation 
regardless of the actual halo mass and its formation history. While this result is encouraging, 
there are some problems. The cooling of gas from the halo should be ongoing, with up to 
lOM yr _1 of cold gas falling onto the disk, well in excess of the observed star formation 
rates (White 1994) and producing X-ray halos much brighter than observed (Fabbiano 1989). 
Even for a minimum steady-state accretion rate of a few M yr -1 , necessary to accumulate a 
typical disk mass in a Hubble time, the dissipation of the radial-infall kinetic energy should 
already produce ~ 10 41 erg s _1 , well beyond the typical limits on diffuse X-ray emission from 
spirals, ^5 X 10 39 ergs _1 (Bregman & Glassgold 1982; Fabbiano 1989). The lack of X-ray 
halos in spiral galaxies suggests that gas accretion from the halo has essentially ended and 
that an isothermal gaseous halo cannot be the regulating mechanism for the Tully-Fisher 
relation. SPH simulations (Navarro & White 1994) support this view by showing that the 
halo gas tends to cool efficiently due to the clumping that characterizes hierarchical merging, 
so that rather little gas is left in X-ray halos at the present time. Although these simulations 
do not include heating from supernova-driven winds that might slow the cooling or return gas 
to the halo, it is unlikely that this heating will maintain the gas just at the virial temperature 
of the halo at all times. 

Empirical evidence that cold gaseous disks have already condensed at early times comes 
from the fact that the mass of neutral hydrogen (HI) in damped Lya systems at z ~ 2 — 3 
is comparable to the total mass in stars of local disk galaxies (Lanzetta et al. 1995). The 
detected HI g column density ~ 10 20 22 cm 2 , similar to that found in nearby galactic 

disks (Broeils & Van Woerden 1994). The damped Lya absorption is often accompanied by 
abundant heavy-element absorption lines at low-ionization stages (Turnshek et al. 1989; 
Wolfe et al. 1993; Petini et al. 1994). The velocity field traced by these metal lines relative 
to the damped Lya line is consistent with typical galactic rotation velocities, as if the damped 
component were associated with a galactic disk while the metal-line systems were embedded 
in a surrounding halo (Lanzetta & Bowen 1992; Lu et al. 1993; Turnshek & Bohlin 1993; 
Lu & Wolfe 1994). The association with galaxies is indeed confirmed by direct imaging in a 
few cases (Steidel et al. 1994, 1995). These facts suggest that a considerable fraction of the 
cold gas necessary to form the stellar content of galactic disks today was already available at 
z ^ 2 and could not have been the result of a recent cooling process. The inferred metalicity 
of this high-redshift gas is an order of magnitude lower than solar (Pei, Fall, & Bechtold 
1991; Sembach et al. 1995), implying subsequent star formation and metal enrichment. If 
the halos of disk galaxies were also formed at z ^ 1 and have not accreted mass afterwards, 
then the low scatter in the Tully-Fisher relation could be explained. 

Evrard, Summers, and Davis (1994) have performed an SPH simulation for an = 1 
CDM model with <r 8 = 0.6. Many of the halos in their simulation had disks of cold gas at 
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the center, and the correlation between the mass of the central disk and the circular velocity 
of the halo proved to have very little scatter, corresponding to ~ 0.13-0.26 mag at a redshift 
of unity. However, since the analysis of the simulated galaxies was performed at z = 1, 
one should use <r 8 = 0.3 in comparing it to the z = spherical model. At the mass scale 
considered, this leads to gm = 1.8 and n = —2, which, according to Figure 1, yields a scatter 
of 0.07 in a or about 0.2 mag in luminosity. Hence, the prediction of the spherical model is 
comparable to that of the simulation. The model also predicts that if the simulation were to 
be evolved down to z = 0, then the scatter in the estimated Tully-Fisher relation would have 
been substantially increased. Future simulations of this type would be particularly valuable 
as they could shed light on some of the issues raised above, including gas cooling, accretion 
in groups of galaxies, and non-sphericity. 

We conclude that the expected scatter in the M-a properties of dark matter halos is 
too large to result in the observed Tully-Fisher relation. We see two paths for breaking the 
assumptions leading to this conclusion. First, if galaxies formed at redshifts z ^ 1 and did 
not accrete material subsequently (contrary to expectations from the standard hierarchical 
clustering picture), then the scatter in the M-a relation would be significantly reduced. 
Second, if a strong feedback process comes into play, it may regularize star formation and 
gas dynamics and detach the overall luminosity of spiral galaxies from details concerning 
the formation history of their halos. Either one of these paths would have significant 
consequences for the process of galaxy formation. The similarity between the Faber- Jackson 
(1976) relation for ellipticals and the Tully-Fisher relation for spirals (cf. Peebles 1993) may 
hint that the actual path taken in reality is universal among different galaxy types. 
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